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■ Abstract 
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m , 

In response to recent experimental studies we investigate Coulomb distortion 



effects on (e, e') reactions from medium and heavy nuclei for the case of elec- 
trons and positrons. We extend our previously reported full DWBA treatment 

O, 

of Coulomb distortions to the case of positrons for the 20S Pb(e, e') reaction 
in the quasielastic region for a particular nuclear model. In addition, we use 
previously reported successful approaches to treating Coulomb corrections in 
an approximate way to calculate the Coulomb distortion effects for (e, e') re- 
actions for both electrons and positrons for the case of a simple nuclear model 
for quasielastic knock-out of nucleons. With these results in hand we develop 
a simple ad-hoc approximation for use in analyzing experiments, and discuss 
methods of extracting the "longitudinal structure function" which enters into 
evaluation of the Coulomb sum rule. These techniques are generally valid 
for lepton induced reactions on nuclei with momentum transfers greater than 
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approximately 300 MeV/c. 
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I. INTRODUCTION 



A persistant problem in using electron scattering for investigating nuclear structure and 
nuclear properties, especially in the quasielastic region, is the large static Coulomb field of 
medium and heavy nuclei. The presence of the static Coulomb potential (of order 25 MeV 
at the surface of the 208 Pb nucleus) invalidates one of the primary attributes of electron 
scattering as usually presented. Namely that in the electron plane-wave Born approximation, 
the cross section can be written as a sum of terms each with a characteristic dependence 
on electron kinematics and containing various bi-linear products of the Fourier transform 
of charge and current matrix elements. That is, various structure functions for the process 
can be extracted from the measured data by so-called Rosenbluth separation methods. The 
trouble with this picture is that when Coulomb distortion of the electron (or positron) 
wavefunctions arising from the static Coulomb field of the target nucleus is included exactly 
by partial wave methods, the structure functions can no longer be extracted from the cross 
section, even in principle. 

In the early 90's, Coulomb distortion for the reactions (e, e') and (e, e'p) in quasielastic 
kinematics was treated exactly by the Ohio University group Jl|-|5| using partial wave ex- 
pansions of the electron wavefunctions. Such partial wave treatments are referred to as the 
distorted wave Born approximation (DWBA) since the static Coulomb distortion is included 
exactly by numerically solving the radial Dirac equation containing the Coulomb potential 
for a finite nuclear charge distribution to obtain the distorted electron wave functions. While 
this calculation permits the comparison of nuclear models to measured cross sections and 
provides an invaluable check on various approximate techniques of including Coulomb dis- 
tortion effects, it is numerically challenging and computation time increases rapidly with 
higher incident electron energy. Furthermore, the initial computer codes did not include the 
option of calculating positron induced reactions in an obvious manner although only the 
sign of the Coulomb distortion term in the Dirac equation needed to be changed. And, as 
noted above, it was not possible to separate the cross section into various terms containing 
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the structure functions and develop insights into the role of various terms in the transition 
charge and current distributions. 

In our DWBA investigations of (e, e') and (e, e'p) reactions in the quasielastic region, 
we used a relativistic treatment based on the a — uj model for the nucleons involved. In 
particular, for the (e, e'p) reaction we use a relativistic Hartree single particle model for a 
bound state |J and a relativistic optical model for an outgoing proton [f?| combined with 
the free space relativistic current operator J M = 7^ + i-^a^ v d v . For the (e, e') case we 
solve for the continuum nucleon wavefunctions using the real bound state potential so as to 
maintain current conservation. Using these models, we compared our DWBA calculations 
with experimental data measured at various laboratories for (e, e') |],|2[|, and for (e, e'p) |3|-[5| 
and have found excellent agreement with the data. We concluded that the relativistic nuclear 
models are in excellent agreement with the measured data and note that we do not need to 
invoke meson exchange effects and other two-body terms in the current that are necessary in 
a Schrodinger description that uses a non-relativistic reduction of the free current operator 
||. However, other investigators use other nuclear models and our elaborate DWBA code 
can not be easily modified to include different transition currents. 

To avoid the numerical difficulties associated with DWBA analyses at higher electron 
energies and to look for a way to still define structure functions, our group [P-[TT[| developed 
an approximate treatment of the Coulomb distortion based on the work of Knoll |12| and the 



work of Lenz and Rosenfelder [13]. We were able to greatly improve some previous attempts 



along this line [|Tj,|Tj] where various additional approximations were made which turned 
out not to be valid. The essence of the approximation is to calculate the four potential 
arising from the lepton four current in the presence of the static Coulomb field of the 
nucleus. This is possible for momentum transfers greater than approximately 300 MeV/c in 
a limited spatial region which we take to be of order 3R where R is the nuclear charge radius. 
The Coulomb distortion is included in the four potential by the elastic scattering lepton 
phase shifts and by letting the magnitude of the lepton momentum include the effect of the 
static Coulomb potential. This last step leads to an r-dependent momentum. A key result 



of our approximation method is that the separation of the cross section into a "longitudinal" 
term and a "transverse" term is still possible. 

We compared our approximate treatment of Coulomb distortion (which we will designate 
as approximate DW) to the exact DWBA results for the reaction (e,e'p) and found good 
agreement (at about the 1-2% level) near the peaks of cross sections even for heavy nuclei 
such as 208 Pb. With an improved parametrization of the elastic scattering electron phase 



shifts ||1 1|| , we achieve quite good agreement away from the peaks in the cross sections. Using 
this approximate DW treatment of Coulomb distortions for the inclusive (e, e') reaction 
in much more difficult numerically since the direction of the outgoing nucleon has to be 
integrated over, and all the nucleons in the nucleus have to be knocked out. Therefore, 
we sought even more severe approximations in order to obtain a simple ad-hoc method of 
calculating the structure functions for (e, e') reactions. In our earlier work, we found it 
necessary to use different ad-hoc procedures for the longitudinal and transverse terms [JTOj] . 



although our investigation of the ad-hoc procedure for the longitudinal terms was hindered by 
the fact that the longtitudinal contributions to the total cross section are usually considerably 
less than 50% and thus we did not have great sensitivity to the Coulomb corrections for 
the longitudinal structure function. In this paper we will use a simple non-relativistic toy 
model to calculate the Coulomb corrections to the longitudinal structure function with 
our approximate DW methods that we applied to (e,e'p) and then investigate the ad-hoc 
treatment of the longitudinal structure function which is a key ingredient in investigating the 
Coulomb sum rule. After developing an improved ad-hoc procedure using our toy model we 
compare it the the full DWBA calculation which we have now extended to include positron 
induced reactions. 

II. APPROXIMATE TREATMENT OF COULOMB DISTORTION 

Our approximate method of including the static Coulomb distortion in the electron 
wavefunctions is to write the wave functions in a plane- wave-like form ||10|; 
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^'( r ) e ±i5(L 2 ) e iA e ip'(r)-r 
P 



(1) 



where the phase factor S(L 2 ) is a function of the square of the orbital angular momentum, 
Up denotes the Dirac spinor, and the local effective momentum p'(r) is given in terms of the 
Coulomb potential of the target nucleus by 



The ad — hoc term A = a[p'(r)-f]L 2 denotes a small higher order correction to the electron 
wave number which we have written in terms of the parameter a = —aZ( 16M ^ v ^ c ) 2 . The 
value of 16 MeV/c was determined by comparison with the exact radial wave functions in 
a partial wave expansion. We have examined the positron case {Z \— » —Z) and find that 
this parametrization works equally well when compared to the exact radial positron wave 
functions. 

We calculate the elastic scattering phases and fit them to a function of the square of 
the Dirac quantum number k used to label the phase shifts. We then replace the discrete 
values of k 2 with the total angular momentum operator J 2 which we subsequently replace 
by the orbital angular momentum operator L 2 since the low k terms where the difference 
between j and I is significant contribute very little to the cross section. Finally we replace 
the angular momentum operator squared by its classical value (r x p) 2 . The removal of any 
spin dependence apart from what is in the Dirac spinor u p is crucial for writing the cross 
section as the sum of a longitudinal and a transverse contribution. 

Initially || we fitted the phases S K to a quadratic function of k 2 which worked reason- 
ably well for lower electron energies, but with the prospect of new higher energy electron 
accelerators, we needed a fit to the phases that will work at higher energies. In addition, 
we wanted to avoid calculating all of the elastic phase shifts, particularly the very high k 
values. We decided to make use of the fact that the higher k phase shifts approach the point 
Coulomb phases which have a simple analytical form at high energy. The low k phases, 
corresponding to orbitals which penetrate the nucleus, are linear in k, 2 which was the basis 




(2) 
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or our initial parameterization. The difficult phases to fit correspond to k values of order 
pR which, from a classical point of view, correspond to scattering from the nuclear surface 
region and are known to make large contributions to the cross section. We were able to 
find a parametrization of the elastic scattering phases shifts in terms of k 2 which has the 
correct large k 2 behaviour and becomes linear in k 2 at low angular momentum, and since 
we have the correct large k behaviour, we need only calculate the exact scattering phase 
shifts for k values up to order pr. After some investigation [fLl|| , we found that the following 
parametrization of elastic scattering phase shift describes the exact phase shifts very well: 

K- 2 i.4>c 2 rv7, k 2 
5(k) = [oo + «2p^]e^ - 2L(1 - e-W)xln(l + k 2 ) (3) 

where p is the electron momentum and we take the nuclear radius to be given by R = 
1.2 A 1 ' 3 — 0.86/A 1 / 3 . We fit the two constants ao and to two of the elastic scattering 
phase shifts (k — 1 and k = Int(pR) +5). To a very good approximation, a = 5(1) and 
a 2 = 4S(Int(pR) + 5) + aZln(2pR), where Int(pR) replaces pR by the integer just less than 
pR. Note that this parametrization only requires the value of the exact scattering phase 
shift for k — 1 and k = Int(pR) + 5. For this paper we have confirmed that this same 
parametrization works equally well for the positron phase shifts. 

Using the new phase shift parametrization and the local effective momentum approxima- 
tion, we construct plane-wave-like wave functions for the incoming and outgoing electrons. 
Since the only spinor dependence is in the Dirac spinor all of the Dirac algebra goes through 
as usual and we end up with a M0ller-like potential given by, 

A appro.DW^ = ^ e _ ^((r X p< (r)?)+5 f ((r xp^ (r)) 2 )] jjAj- A f ) c »q'(r)-r ^ ^ 

where the phase shift parametrization is given in Eq. 3 with k 2 being replaced by (r x p) 2 , 
the parameter A is given following Eq. ^|, and the r-dependent momentum transfer is given 
by q'(r) = pj(r) - P/(r). 

With this approximate DW four potential A^ it is straightforward to calculate the (e, e'p) 



cross sections and modified structure functions. We showed [jlT] that using this new phase 



shift(see @,[T(J for details) we can reproduce the full DWBA cross sections for (e, e'p) from 
medium and heavy nuclei very well. 



III. APPLICATION TO THE INCLUSIVE PROCESS 

In the plane wave Born approximation (PWBA), where electrons or positrons are de- 
scribed as Dirac plane waves, the cross section for inclusive quasielastic (e, e') processes can 
be written simply as 

i]!(J a M { q 4s L (q,w) + [tan 2 % - f,]S T (q,w)} (5) 



where g 2 = uo 2 — q 2 is the four-momentum transfer, Om is the Mott cross section given 
by cr M = (^) 2 COS 4 1 , and Sl and St are the longitudinal and transverse structure functions 

sm 2 

which depend only on the momentum transfer q and the energy transfer u. As is well known, 
by keeping the momentum and energy transfers fixed while varying the electron energy E 
and scattering angle 6 e , it is possible to extract the two structure functions with two mea- 
surements. As we will summarize below, our approximate treatment of Coulomb distortions 
still permit Rosenbluth-like separations but with Coulomb corrections which require the use 
of models. 

For the inclusive cross section (e, e'), the longitudinal and transverse structure functions 
in Eq. ([5]) are bi-linear products of the Fourier transform of the components of the nuclear 
transition current density integrated over outgoing nucleon angles. Explicitly, the structure 
functions for knocking out nucleons from a shell with angular momentum jj, are given by 



StM = £ / (|Ag 2 + \N s ?)d{l„ (7) 

where the nucleon density of states p p = the z-axis is taken to be along q, and 

and s p are the z-components of the angular momentum of the bound and continuum state 
particles. The Fourier transfer of the nuclear current J^(r) is simply, 
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= / J^(r)e iqr rf 3 r. (8) 



ans the continuity equation has been used to eliminate the ^-component (N z ) via the equa- 



tion N z = -ajv . 



When we use our approximate M0ller potential given in Eq. (|]) , we also can separate 
the cross section into longitudinal and tranverse components since as noted previously, it 
the Dirac spinor structure that leads to this result. However, when we use the approximate 
electron four potential along with current conservation to eliminate the z-component of 
the current we run into a problem since the momentum transfer q' depends on r both in 
magnitude and direction. In addition, the phase factors depend on r. To avoid generating 
additional terms we assume the direction of q'(r) is along the asymptotic momentum transfer 
q which defines the z-axis, and neglect the dependence on r in the phases and in q'(r), when 
taking the divergence of N. With this further approximation, current conservation implies 
lvNq + q'(r)-N = 0. Using these results, approximate the cross section for the inclusive 
reaction (e, e') can be written as 

_g_ = <TM{ | 5l( ,, m ) + [tan2 |_^ ]5 ; ( ,, m)} (9 ) 

and the transform of the transition nuclear current elements which appears in Sl and St 
are given by 



N 



appro. DW 



(<M.f( ? )2 e «/([rxpJ(r)]») e tfA[«^(r)n e i(A i -A / ) e iq'(r).r Jo(r)d 3 r (1Q) 

% q r (r) 

■^appro.DW _ J ^([rxp'^r)] 2 ) e i5 f ([rxp' f (r)} 2 ) ^(A.-Af) e iq_' (r)-r J T ^£ r (11) 

Due to the angular dependence in the phase factors in Eqs. (|TUj) and (|TT|), a multipole 
expansion of the approximate potential is not practical. Thus, jy^ppro.DW an( ^ -^appro.DW 
to be evaluated by carrying out a 3-dimensional numerical integration. As we have shown 
for the (e, e'p) case [EC], this numerical integration reproduces the exact DWBA results very 
well. However, since the inclusive reaction (e, e') requires a sum over all occupied neutron and 
proton shells and a further integration over the directions of the outgoing nucleon, numerical 
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integration is very time consuming. In order to have a more practical procedure we examine 
additional approximations that will allow the integration over the angular coordinates in 



Eqs. ( |I0D and (|TT| ) to be done analytically. 

We created such an ad-hoc procedure in a previous paper ||, but we were comparing 
our ad-hoc procedures to the exact DWBA calculation which was largely dominated by the 
transverse terms. Hence, our ad-hoc procedures for the longitudinal term were not very well 
determined. In addition, our full DWBA calculation was only set up for electrons, so we 
could not check the ad-hoc approximation for positrons. In order to address this matter, 
we created a simple toy model which assumes harmonic oscillator bound state protons and 
takes the outgoing continuum proton wavefunction to be a plane wave. Using this simple 
model to calculate the transition charge distributions allows us to calculate the longitudinal 
contribution to the cross section using the approximate DW expression for N Q of Eq. (|10"D 
and to compare this result to various ad-hoc proscriptions. Based on this investigation, 
coupled with our previous investigation of the transverse contributions which dominate the 
cross section at large electron scattering angles, we propose the following ad-hoc expressions 
for the longitudinal and transverse structure functions: 



nrad—hoc 
JV 



-^ad—hoc 

T Pi ' 



,^0)2, Q )2 e i<«(« i a )+^)> e iq'(r).r Jo(r)d 3 r ^ 

J e iq ' (r) - r J T (r)d 3 r. (13) 



where < 5(k 2 j) > denotes an average over the angles of the vector r. That is, < k\ j >=< 
(r x Pi,/) 2 >= r 2 pfj(3 — cos 2 9 Pi f )/4. Note that under this averaging, the A term goes 
to zero. This removes the angular dependence in the phase factors, and thus permits a 
multipole treatment of the matrix element as usual. 

In the following figures for the longitudinal parts of the cross sections based on our simple 
model we will compare our new recommended longitudinal ad-hoc result given in Eq. fll2|) 
to the result calculated by the full three-dimensional integration of Eq. ( fLOf) and to our 
previous ad-hoc results called LEMA' which we give below for convenience: 
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N LEMA> = ( ^(0) f e V'(r).r J()(r)rf 3 r (M) 



where q" = p"(r) ~~ P/( r ); £>"( r ) = V ~~ 7 Jo ^( r ')^ r ' an d the factor A, which depends on the 
energy transfer u, is given by A = (uj/uj ) 2 with u = r|rr- 

Clearly Nq d ~ hoc and N^ d_/loc represent a modified Fourier transform of the nuclear tran- 
sition current. For comparison purposes, the approximation known as the EMA replaces 
q'(r) with q'(0) wherever it appears in Eqs. ( |T0"D and ( II ) for 7V"o and Nt and the phases 
are neglected as usual. We find that for light nuclei the EMA is adequate, but it leads to 
large errors for nuclei as heavy as 208 Pb. 

In Fig. 1 we compare the two approximate calculations with the DW approximation for 
the longitudinal contribution to the cross section for knocking protons out of various shells 
at a forward angle in 208 Pb by electrons or positrons. Note that while we use harmonic 
oscillator wavefunctions for all orbitals, we do use the binding energies of the orbitals that 
correspond to the values we find for our relativistic a — uj model for 208 Pb. While the ad-hoc 
result is not in perfect agreement with the full DW result, it clearly is in better agreement 
that the LEMA' result and, for cases where the electron incident and final energy exceed 
300 MeV is in reasonable agreement, particularly near the maxima. Note that the positron 
results are not very sensitive to which approximation is used. 

In Fig. 2 we show similar results at a backward angle. We note that our ad-hoc DWBA 
results for positrons tend to be in much better agreement with the DW result than the 
electron case. We again find that our new ad-hoc approximation for the longitudinal con- 
tribution is considerably better that our previous LEMA' result. We note that while the 
agreement between our ad-hoc calculation and the DW calculation for knocking out protons 
from individual orbitals is not excellent, the discrepancies do not seem have a systematic 
tendency to be either low or high and we have reason to hope that when all the orbitals are 
added together as in the case of (e, e') reactions from nuclei that these discrepancies will 
tend to average out. 



11 



IV. COMPARISON TO EXPERIMENT AND CONCLUSIONS 



Based of our investigation of this simple toy model, we adopt our new ad-hoc model for 
the longitudinal structure functions and return to our full nuclear model for investigating 
Coulomb corrections for 208 Pb(e, e') in the quasielastic region where the lepton can be elec- 
trons or positrons. Our first step is to re-examine our full DWBA calculation || and modify 
the code for the case of positrons. We were successful in doing this and can now compare the 
full DWBA calculation for electrons and positrons based on a realistic relativisitic nuclear 
model to our ad-hoc treatment of Coulomb corrections which still permit a separation into 
longitudinal and transverse terms. 

In Fig. 3 we compare the full DWBA calculation to the ad-hoc result and the electron 
plane wave result (PWBA) for electrons with incident energy of — 310 MeV and scattering 
angle of 9 = 143°. Note that this comparison is only a test of the our ad-hoc transverse 
treatment which is unchanged from our previous work since the longitudinal contribution at 
such a large angle is at the few percent level. The agreement of the " plane- wave-like" ad-hoc 
calculation with the full DWBA result is quite good, even though the outgoing electron 
energy is well below 300 MeV. 

In Fig. 4, we perform a similar comparison for the case of positrons with incident energy 
of Ei = 485 MeV and scattering angle 9 = 60°. Again, the agreement is quite good, and 
unlike the backward angle scattering case, the longitudinal response contributes about 40% 
of the cross section. 

We have examined a number of other cases, and the agreement shown in Figs. 3 and 4 
is characteristic at these energies. As the lepton energies increase, the ad-hoc approxima- 
tion improves since the Coulomb distortion effects become smaller. We did notice in our 
investigations a general tendency that Coulomb distortion effects for positrons tend to be 
smaller than Coulomb distortion effects for electrons. This corresponds to an observation 



made many years ago when looking at inelastic lepton scattering from nuclei [|16|| , where 
we noted that Coulomb distortion for positrons tends to saturate. As electrons pass near 
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the nucleus, the attactive Coulomb interaction pulls them into regions of stronger potential 
which increases the Coulomb distortion effects, while positrons are pushed away from the 
region with a stronger potential. 

With our capability of examining Coulomb distortion of both positrons and leptons 
with the full DWBA calculation and with our improved ad-hoc procedure we can compare 
our model predictions to experiment. In Fig. 5, we compare our model calculations with 
Coulomb distortion included exactly and with our ad-hoc method for quasielastic scattering 
of electrons of energy 383 MeV and positrons of energy 420 MeV both at a scattering angle 
of 9 = 60° from 208 Pb to the experimental data from Saclay [0,0. Note that in this and 



the following figure, we are plotting the total structure function S to tai = ^fm /a M {Ei). 



We first note that our ad-hoc and exact DWBA results are in reasonable agreement 
although the lepton energy is somewhat low for our approximate result, and further that 
the positron and electron total structure functions have approximately the same shape as 
a function of the energy transfer u. However, they do not have the same magnitude as 
do the data from Saclay. The positron theory result is in reasonable agreement with the 
experimental data, but the electron result is approximately 15%-20% larger than the data. 

In Fig. 6 we make a similar comparison except that now the scattering angle is 6 = 
143°, and the electron incident energy is 224 MeV while the positron incident energy is 
262 MeV. Again, when Stotai is plotted the positron and electron shapes as a function of 
energy transfer u are very similar, but again, unlike the experimental data, the magnitudes 
are quite different. At this backward scattering angle case, our electron result (DWBA) is 
in quite good agreement with the data. At these much lower energies, clearly our ad-hoc 
approximation is beginning to fail, particularly for the electron case. 

There is considerable interest in extracting the longitudinal contributions from (e, e') re- 
actions from medium and heavy nuclei in order to investigate the Coulomb sum rule. Clearly, 
Coulomb distortion effects have to be handled properly. Our results indicate that we could 
use a Rosenbluth-like procedure in order to separate our "longitudinal" and "transverse" 
contributions to the cross section. However, these contributions depend on a modified (by 

13 



Coulomb distortion) Fourier transform of the transition charge and current distributions. It 
is necessary to use a nuclear model to extract the longitudinal and/or tranverse structure 
functions from the data. It is not clear to us that a Rosenbluth-like procedure is the best 
way to proceed, since our ad-hoc procedure is not accurate in the wings of the cross section 
distributions and in many cases, some of the Rosenbluth points fall on either the low uj or 
high uj side of the quasielastic peak. It seems that a better procedure might be to choose 
some semi-realistic nuclear model for the process in question. Use Eqs. (|l^) and (|13|) to 



calculate the structure functions and then fit the calculations to the available data using a 
least squares procedure to determine nomalization factors Nl and Nt in front of the ap- 
propriate terms. The nuclear model should have the overall correct spatial and kinematic 
dependence, but the longitudinal or transverse strength will be determined by fitting these 
normalization factors. Once these factors are determined, one can use the same nuclear 
model weighted with these factors to calcuate the plane wave structure functions, thereby 
having "measured" the nuclear longitudinal and transverse response. 
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FIGURES 

FIG. 1. Longitudinal contributions to the differential cross sections at a forward scattering angle 
for 20S Pb(e±, e±) for different bound state orbitals. The solid line is the approximate DW result, 
the dashed line is our ad-hoc result and the dotted line is our previous LEMA' approximation. 

FIG. 2. Longitudinal contributions to the differential cross sections for 208 P6(e-i-, e^_) at a back- 
ward angle for different bound state orbitals. The solid line is the approximate DW result, the 
dashed line is our ad-hoc DWBA result and the dotted line is our previous LEMA' approximation. 

FIG. 3. The DWBA differential cross section for 208 P6(e_, e'_) at 310 MeV and scattering 
angle 9 = 143° compared to our ad-hoc DWBA and to the plane wave result. The bound state and 
continuum neutron and proton orbitals are solutions to relativistic Hartree potential based on the 
a — uj model. 

FIG. 4. The DWBA differential cross section for 208 Pb(e + , e' + ) at 485 MeV and scattering 
angle 9 = 60° compared to our ad-hoc DWBA and to the plane wave result. The bound state and 
continuum neutron and proton orbitals are solutions to relativistic Hartree potential based on the 
a — uj model. 

FIG. 5. The total structure function S 'total generated by dividing the differential cross section 
by o~m for 208 P6(e±, e'±) at a forward scattering angle of 60° with electrons of energy 383 MeV and 
positrons with energy 420 MeV. The theoretical curves correspond to the full DWBA calculation 
and to our ad-hoc DWBA calculation. The data were taken at Saclay [17,18|. The bound state 



and continuum neutron and proton orbitals are solutions to relativistic Hartree potential based on 
the a — uj model. 
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FIG. 6. The total structure function Stated generated by dividing the differential cross section by 
om for 20S Pb(e±, e'±) at a backward scattering angle of 143° with electrons of energy 224 MeV and 
positrons with energy 262 MeV. The theoretical curves correspond to the full DWBA calculation 



and to our ad-hoc DWBA calculation. The data were taken at Saclay [17,18|. The bound state 
and continuum neutron and proton orbitals are solutions to relativistic Hartree potential based on 
the a — to model. 
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